Non-zero transversity distribution of the pion in a quark-spectator-antiquark model 
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We calculate the non-zero (naive) T-odd transverse momentum dependent transversity distribu- 
tion /ii~(a;,kj_) of the pion in a quark-spectator-antiquark model. The final-state interaction is mod- 
elled by the approximation of one gluon exchange between the quark and the antiquark spectator. 
Using our model result we estimate the unsuppressed cos2<^) azimuthal asymmetry in unpolarized 
ir~p Drell-Yan process. We find that the transverse momentum dependence of h~t (x, kjj of the 
pion is the same as that of h\{x, k^_) of the proton calculated from the quark-scalar-diquark model, 
although the x dependencies of them are different from each other. This suggests a connection 
between cos2<^) asymmetries in Drell-Yan processes with different initial hadrons. 
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I. INTRODUCTION 

Recently the study of transverse momentum depen- 
dent distribution functions is among the special is- 
sues in hadronic physics. Of particular interest, are 
two leading-twist (naive) time-reversal odd transverse 
momentum dependent distribution functions: Sivers 
function [ll Q fyr(x, k^_) and its chiral-odd partner 
/ii~(x,kj_) 0, 3- Sivers function represents the unpo- 
larized parton distribution in a transversely polarized 
hadron, while (x, kj_) denotes the parton transversity 
distribution in an unpolarized hadron. One main moti- 
vation to investigate these two distributions is that they 
are the possible sources of the unsuppressed azimuthal 
asymmetries observed in hadronic reactions. The former 
distribution function was proposed first by Sivers |l( to 
illustrate that it can lead to large single-spin azimuthal 
asymmetries. This nontrivial correlation between the 
transverse momentum of the quark and the polariza- 
tion of the hadron was thought to be forbidden by time- 
reversal invariance Recently a direct calculation 0,0 
of Sivers asymmetry by inclusion of final-state (in semi- 
inclusive deep inelastic scattering(SIDIS)) or initial-state 
interaction (in Drell-Yan process) shows that the asym- 
metry is in principle non-zero. Then it was found that 
the presence of the Wilson lines in the operators defining 
the parton densities allows for the Sivers effect without 
a violation of time-reversal invariance |8|, and the final- 
or initial-state interaction can be factorized into a full 
gauge-invariance definition of transverse momentum de- 
pendent distribution functions 0, [^} . 

These theoretical developments open a wide range of 
phenomenolo gica l applications. Several model calcula- 
tions 0, H2l Il3l ITl ] of Sivers function have been per- 



formed to estimate single-spin asymmetries in SIDIS 
process, which is under investigation by current exper- 
iment ^5(. On the other hand, it is shown that 
non-zero h^(x, kj_) can arise from the same mecha- 
nism which produces /j^(x,k^). It has been demon- 
strated Q| that hi(x,k\) can account for the substan- 
tial cos2(/) asymmetries in unpolarized Drell-Yan lepton 
pair production from pion-nucleon scattering: N — > 
H + /i~X [H,II3- In Ref. [3, Boer, Brodsky, and Hwang 
computed h^(x, kj_) of the proton in a quark-scalar- 
diquark model within soft gluon exchange. They found 
that hi (a;,kj_) is equal to f{^(x, kj_) obtained from the 
same model. Then the maximum magnitude of the cos2</> 
asymmetries mpp — > UX is estimated to be ^30%, by us- 
ing the calculated hj-(x,\i.\). In this paper, we perform 
the first computation on hj~(x, kj_) of the pion (denoted 
as h^(x, k^_)) in a quark-spectator-antiquark model in 
presence of final-state interaction, in similar to the quark- 
scalar-diquark model of the proton. We find that the 
transverse momentum dependence of hi 7T (x,\i 2 L ) in our 
model is the same as that of h^{x,\i 2 1 ) of the proton 
from the quark-scalar-diquark model, although the x de- 
pendence is different. This feature allows one to expect 
that hi v (x,'k 2 L ) and /i;f (x,k^_) are closely related. With 
the present model result we investigate the cos20 asym- 
metry in unpolarized n~p Drell-Yan process and obtain 
an unsuppressed result. The shape of the asymmetry is 
similar to the cos20 asymmetries in pp — > UX, estimated 
in Ref. 01- The result suggests a connection between 
cos20 asymmetries in Drell-Yan processes with different 
initial hadrons. 



II. CALCULATION OF ft^(a;,ki) OF THE PION 



'Corresponding author. 
mabq@phy.pku.edu.cn 



Electronic address: 



In this section, we present the calculation /^(a^k^) 
of the pion. We start our computation from the quark 
light-cone correlation function of the pion in Feynman 
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gauge (we will perform our calculation in this gauge): 



$ Q/3 (a;,kj 



(2tt) 3 



P ff 1^(0)4(0", oo-) 



(1) 



We use notation a ± = a° ± a 3 , a ■ b — ^(a + b + 
a~b + ) — • bj^. The pion momentum is denoted by 
f£ = (P+,p-,P n± ) = (P+,M/P+,0 X ). £ (0,oo) is 
the path-ordered exponential (Wilson line) with the form: 
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FIG. 1: Diagram which gives <£> of the pion in the quark- 
spectator-antiquark model in lowest order. 



£ (0, oo) = V exp -ig / dC • A(0, £ 



(2) 



Releasing the constraint of (naive) time-reversal invari- 
ance and keeping parity invariance and hermiticity, the 
quark correlation function of the pion can be parameter- 
ized into a set of transverse momentum dependent dis- 
tribution functions in leading twist as follows 0, fli^ 



<&(cc,kj 



/i(*,ki) a - k ^ 



( 3 ) 

where n is the light-like vector with components 
(n+n-, nj.) = (1,0, Ox), a„ v = §[7",7l, fm(x,kl) 
and hi v (x, kj_) denote the unpolarized quark distribu- 
tion and the quark transversity of the pion, respectively. 
Knowing <$>. K (x,'k.±), one can obtain these distributions 
from equations: 



= Tr[$(a;,k ± )7+]; 
= Tr[$(a:,k L )o i +] 



(4) 
(5) 



We will calculate above distribution functions in the 
quark-spectator-antiquark model. It is similar to the 
quark-scalar-diquark model for calculating Sivers func- 
tion and of the proton, and the differences are that 
the intermediate state here is the constituent antiquark 
instead, as shown in Fig. 1, and that the pion-quark- 
antiquark interaction is modelled by pseudoscalar cou- 
pling: 



(6) 



in which g n is the pion-quark-antiquark coupling con- 
stant, and e2 is the charge of the antiquark. fi w (x,kj_) 
can be calculated from the lowest order <J>(x, kj_) without 
the path-ordered exponential. From Fig.^we obtain 
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where m is the mass of the outgoing quark which is 
the same of the antiquark, v s is the spinor of the an- 
tiquark and the quark momentum k — {xP + , (kj_ + 



m 2 )/xP + , kj_). We take the spin sum as ^ s v s v s = 
(Ifn — ft — m), which is a little different from the spin 
sum adopted in Ref . [2(j . Immediately we arrive at 



/i^(x,ki) = 



{1 - x)[kl + (1 + x) 2 



2(2tt) 3 (k 2 ± + m 2 - x(l - x)MkY 



— 



(ki+B.) 2 ' 



(8) 



where we set CV = - x)/[2(2tt) 3 ], D„ = (1 + x) 2 m 2 
and 



B n = m 2 -x(l-x)Ml. 



(9) 



The T-odd distribution hi w (x, kj_), however, is absent 
in the lowest order kj_). In order to produce this 
T-odd distribution, the path-ordered exponential which 
ensures gauge invariance of the distribution function has 
to be included. The exponential serves as the final-state 
interaction (FSI) or initial-state interaction (ISI) between 
the struck quark and the remnant of the hadron, which 
is also viewed as the soft gluon scattering, to provide 
nontrivial phase to generate T-odd distribution function. 
In our calculation we expand path-ordered exponential 
to hrst order, means that the final- or/and initial-state 
interaction is modelled by one gluon exchange, as shown 
in Fig. Thus the nonzero fij^fe, kj_) can be calculated 
from the expression 



• i -><■■! ^ A+(o,c,o±)dc j oyMOW*)^ 

+h.c, (10) 



in which \q) represents the antiquark spectator state, and 
e± is the charge of the struck quark. In momentum space 
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FIG. 2: Diagram which yields $ with final-state interaction 
modelled by one gluon exchange. 



we write down: 



2h^(x,k 2 L )k i L _ -ieie-2 \p f d 4 q _ t 
M n ~ 8(2tt) 3 (1 - x)P+^J (2tt)* V 

k 2 — m 2 (k + q) 2 — m 2 ^ 5 (k + q + P^) 2 + m 2 
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h.c. 



(11) 



The 7 + in the second line of Eq. comes from the 
antiquark-gluon interaction vertex. The q~ integral can 
be done from the contour method. q + integral is realized 
from taking the imaginary part of the eikonal propagator: 
1/(<Z + + ie) — > — iir8(q + ). Then we obtain 



5^eie 2 (l -x) f d 2 q± 



2(2^(^+^)7 (2tt) 2 

2mq l ± 
: 9l[(k± + q±) 2 + B„Y 



(12) 



To arrive at above equation, we have calculated the 
trace in the numerator in Eq. (|llfl as follows 



+„,« 



(13) 



= Tr[(# w - t - m)(# - m) 7 * 7 +(# + 

x(^ + ^-^ + m) 7 +] (14) 
= 8(P+) 2 rnqY when <?+ = 0. (15) 

After integrating over q^, we obtain the expression of 
hin(x, kj_) in the antiquark spectator model: 



gl eie 2 mM 7r (l - x) fk^ + S w 



2(2^) 3 4tt k^ki + B^) V 



-In 



with iJ^ given in Eq. © and 



(16) 
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In Ref. 0| /^(x, kj_) of the proton is computed in the 
quark-scalar-diquark model as: 



hi(x,k 2 



K(K + b p ) 



in 



B r 



B„ 



(18) 



Comparing Eq. 1)16(1 with Eq. (|18[) one can find that 
the calculated hj^(x, kj_) in the spectator antiquark 
model has a same transverse momentum dependence of 

(x, k^_) of the proton obtained from the quark-scalar- 
diquark model, although the x dependence is different in 
A p /n and B p / V respectively, due to the different mass pa- 
rameters in them. This feature may not be held generally, 
but one can expect that there is close relation between 
distributions of different hadrons since both functions 
are generated by the same underling mechanism. 

We can also calculate the T-odd distribution of the 
transversity distribution of the valence antiquark hj^ of 
the pion from the same model, by replacing the antiquark 
spectator with the quark spectator. A similar calculation 
yields hj^ — h^, showing that the magnitudes of the 
distributions for two valence quarks are the same. 

With fin(x, kjj and hj^{x, k\) we estimate 
\k±hi„(x,\t\)\/[M„fin(x,\s. 2 L )], which is proportional to 
the cos2</> asymmetries in unpolarized Drell-Yan process. 
We choose the mass parameters m = 0.1 GeV for the 
constituent quark mass and M r = 0.137 GeV for the 
pion mass. To generalize our model result to the con- 
sequence in QCD we extrapolate the coupling constant 
eie2/47r — > Cfol s , and take a s = 0.3 which is adopted 
in Ref. 0. We plot the x dependence of the ratio at 
|kj_| = 0.3 GeV in Fig. and the |kj_| dependence 
at x = 0.15 in Fig. |3Jd. The ratios are comparable 
in magnitude with \k±hi(x 7 k 2 ± )\/[M fi(x, k 2 ± )] of the 
proton (one can see Ref. Q, where |kj_/j^|/(M/i) is 
given, since hi = fix) in the quark-scalar-diquark 
model, where M is the proton mass and /i(x,kj_) is the 
unpolarized quark distribution of the proton. 



III. COS2<^ ASYMMETRIES IN UNPOLARIZED 
DRELL-YAN PROCESS 

The general form of the angular differential cross sec- 
tion for unpolarized n~p Drell-Yan process is 



1 da 



1 



4tt A + 3 



1 + Acos 2 6> + /isin #cos0 



-—sin (9cos20J , 



(19) 



where <fi is the angle between the lepton plane and the 
plane of the incident hadrons in the lepton pair center 
of mass frame (see Fig.0J). The experimental data show 
large value of v near to 30%, which can not be explained 
by perturbative QCD. Many theoretical approaches have 
been proposed to interpret the experimental data, such 
as high- twist effect [2l], 123 . and factorization breaking 
mechanism [23] . In Ref. 4] Boer demonstrated that un- 
suppressed cos20 asymmetries can arise from a product 
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FIG. 3: Model prediction of |kj_ft^(a;,ki)|/[M w /i w (x,kj_)] as a function of x and |kj_|. 



of two chiral-odd hj~ which depends on transverse mo- 
mentum. In Ref. |18| the cos2</> asymmetry in unpo- 
larized pp — > HX Drell-Yan process has been estimated 
from h±(x, kj_) for the proton computed by quark-scalar- 
diquark model. The maximum of v in that case is in the 
order of 30%. 

In this section we give a simple estimate of cos20 asym- 
metry in unpolarized 7r~p Drell-Yan process, from ft^ 
computed by our model. The leading order unpolarized 
Drell-Yan cross section expressed in the Collins-Soper 
frame |2J] is 



km 



xcos20JF 



(2h • Px h • k x ) - (p± • k ± ) 



M 1 M 2 



(20) 



where Q 2 = q 2 is the invariance mass square of the lepton 
pair, and the vector h = qx/Qr- We have used the 
notation 



d 2 P± d 2 k ± 6 2 (p ± +kx - qx)/?(z,Px) 
x/f(2,ki). (21) 



From Ea. 1201 one can give the expression for the asym- 
metry coefficient v 



a, a 
a , cc 

_ 2 Ea,a e a^a 

"~ M,M 2 j: ara e 2 G a - 



(2h • pxh • kx) - (px • kj 



M 1 M 2 



(22) 




FIG. 4: Angular definitions of unpolarized Drell-Yan process 
in the lepton pair center of mass frame. 



Our model calculation has shown = hj^.. Thus 
in n~p unpolarized Drell-Yan process we can assume it- 
quark dominance, which means the main contribution to 
asymmetry comes from E l7 j"(^,kj_) x h± ' u (x, Pj_)> since 
u in 7r~ and u in proton are both valence quarks. Then 
we have v w 2F U /(M„MG U ). To evaluate i/, we use our 
model result for h^ u and fi^, and we adopt h^' u and 
/i from Ref. Using the px integration to eliminate 
the delta function in the denominator and numerator in 
Eq. (122H one arrives at 



F u = / d|kx| / d0 qk 
Jo Jo 



! ki(k»+iV 



(k ± - 2k ± cos% fc + |kx||qx|cos0 9fc ) 



xln 



(k 2 ± +/)(ki + / + i3 p ) 
ki + / + B p 



B 



G u = [ d\k±\ f 
Jo Jo 



p 

2n 



d9, 



qk ~ 



q c p |kx|(kj + ^) 
(ki + B^) 2 



(kj + / + D p ) 
(k 2 ± +f + B p ) 2 ' 



(23) 



(24) 



where / = q ± — 2|qx||kx|cos6', jfe , and 8 qk is the angle 
between kx and qx- In above equations we change the 
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FIG. 5: Numerical result for the cos2<^) asymmetry in iv~p 
unpolarized Drell-Yan (solid line), using Al^ — 0.137 GeV, 
m=0.1 GeV, and x — x — 0.2. The dashed line corresponds 
to the same asymmetry in pp Drell-Yan process estimated in 
the quark-scalar-diquark model (see Ref. 01)- 



integration variables into polar coordinates. For the pa- 
rameters in Eq. [33] and Eq. [2] we take A p /C p = 0.3 
(GeV 2 ) and D p = AB p w 1/4 which are used in Ref. [H|, 
and Ax/Ci; = 0.01 (GeV 2 ), = 2B n « 0.014 which 
corresponds x — 0.2. We give the numerical evaluation 
of the asymmetry as a function of Qt in Fig. [3] which 
shows an unsuppressed result of the order of O (10%). 
Besides this, the shape of the asymmetry is similar to 
the cos20 asymmetry in pp unpolarized Drell-Yan pro- 
cess estimated in the quark-scalar-diquark model |l8j . 



the pion for the first time in a quark-spectator-antiquark 
model in presence of final-state interaction, and investi- 
gate the cos2</> asymmetry in unpolarized Drell-Yan pro- 
cess with our model result. The calculated h^(x, kj_) 
shows a same form of transverse momentum depen- 
dence as that of h^{x, kj_) for the proton computed 
from a quark-scalar-diquark model. The similarity of 
h^{x, kjj for different hadrons (for example, nucleon 
and pseudoscalar meson) implies that these functions 
are closely related, because the mechanism that gen- 
erates them is the same. Besides this, hi„(x, kj_) is 
an interesting observable that can account for the large 
cos20 asymmetry in n~ N unpolarized Drell-Yan pro- 
cess. The contributed asymmetry of /i^(x, kj_) is pro- 
portional to \k±hi„\/(M 7r fi 7r ), which is comparable with 
|kx/ii |/(M/i) of the proton in magnitude. With nonzero 
hi v (x, kj_) we reveal an unsuppressed cos20 asymme- 
try in unpolarized n~p Drell-Yan process from our cal- 
culation. The shape of the asymmetry is similar to 
that of cos20 asymmetries in pp Drell-Yan process esti- 
mated in Ref. [lg , suggesting a connection between cos2<^> 
asymmetries in Drell-Yan processes with different initial 
hadrons. 
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